Velocity Inversion in Nanochannel Flow 
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The nanoscale cylindrical Couette flow is investigated by means of molecular dynamics simula- 
tions, in the case where the inner cylinder is rotating whereas the outer cylinder is at rest. We find 
that the tangential velocity of the flow is inverted when the fluid-wall interaction near the outer 
cylinder is weak and the fluid density is low. The unusual velocity inversion behavior is shown 
to be strongly related to the degree of the slip between the fluid and the outer cylinder, which is 
determined by the presence or absence of the layering of the fluid near the outer wall. 

PACS numbers: 61.20.Ja, 83.50.Ha, 83.50.Lh, 83.50.Rp 
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I. INTRODUCTION 



Couette flow between two concentric rotating cylinders 
shows interesting and unexpected feature such as "veloc- 
ity inversion" which implies that the tangential velocity 
of the flow increases with distance from a rotating cylin- 
der to a stationary cylinder. The velocity inversion phe- 
nomenon has been studied with analytical and numeri- 
cal methods. Einzel et al.^ first predicted by suggest- 
ing a generalized slip boundary condition for incompress- 
ible flow over curved or rough surfaces that the velocity 
profile would become inverted in the case of large veloc- 
ity slip at the wall surfaces. Tibbs et al.fy and Aoki 
et o/.[3( confirmed the prediction of the anomalous be- 
havior using direct simulation Monte Carlo calculations 
and some analytic approaches and found that the values 
of the accommodation coefficients in their models play 
an important role on the velocity inversion phenomenon. 
In all those studies, they found that the velocity inver- 
sion phenomenon only occurred for large velocity slip at 
small value of the accommodation coefficient. The ac- 
commodation coefficient is an important parameter in 
determining the degree of slip of the fluid at the wall and 
represents the average tangential momentum exchange 
between the flowing fluid particles and the wall bound- 
ary. The case of zero value of the coefficient is called 
specular reflection, meaning zero friction. When its value 
is unity, the reflection is diffuse, meaning that the fluid 
particles are reflected with zero average tangential veloc- 
ity. Recently, similar studies for the velocity inversion 
phenomenon were performed [H, [H, Q . It, however, is dif- 
ficult that the concept of an accommodation coefficient 
can apply to all types of fluids or the nature of the wall 
materials. 

Molecular Dynamics (MD) is a very powerful tool in 
the exploration and study of the nature of the flow at 
the boundaries independent of the properties of the fluid 
and wall materials. For simple Lennard-Jones liquids, 
the wetting properties of the fluid can be modeled by 
wetting parameter varying the strength of the fluid-wall 
attraction [7], la |9(. The wetting parameter could directly 
be related to the accommodation coefficient of the fluid- 
wall system. When the value of parameter is zero (non- 



wetting or specular reflective), the fluid- wall interface is 
specular reflective similar to the case of zero value of the 
accommodation coefficient. When it is unity (wetting or 
attractive), the surface is diffuse reflection. The large 
slip at the fluid-wall interface exhibits under nonwetting 
condition and it decreases as the value of wetting pa- 
rameter increases to unity. The degree of the slip is also 
dependant on the density of the fluid. Even at a strongly 
attractive wall, a large slip was found in the highly dilute 
gas regime 0. 

In this paper we use MD simulations for Lennard- Jones 
liquids to investigate the velocity inversion phenomenon 
in nanoscale cylindrical Couette flows of the concentric 
rotating cylinders, varying the fluid density and the value 
of the wetting parameter. The tangential velocity of the 
flow is measured in the case where the inner cylinder is 
rotating whereas the outer cylinder is at rest. At low 
fluid density and under the poor wetting condition, the 
velocity is inverted with large velocity slip near the outer 
cylinder. On the other hand, at high density or under 
the good wetting condition, such inversion does not oc- 
cur with the layer adsorption of the fluid near the outer 
cylinder. Next, the radial force field is considered, and 
it is found that the radial force is a better measurement 
than the fluid density to examine the effect of the layering 
near the outer cylinder on the velocity slip. 




FIG. 1: (Color online) Views of a fluid in a nanochannel con- 
sisting of two concentric cylinders. 
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II. MOLECULAR DYNAMICS SIMULATIONS 

To investigate the properties of the nanochannel flows 
in a curved surface, we have conducted MD simulations 
of the cylindrical Couette flows in the concentric rotat- 
ing cylinders as shown in Fig. [T] Standard MD tech- 
niques are used [l(| E|- The fluid particles and the wall 
particles at a distance r interact through a Lennard- 



Jones(LJ) potential, Vij(r) — 4e[(r/cr) 



-12 



A tJ (r/aY 



where e and a represent the energy and length scales, 
and Aij = Aji is a dimensionless parameter that con- 
trols the attractive part of the potential for the fluid-fluid 
and the fluid-wall particles. As mentioned previously, 
this parameter is similar to the accommodation coeffi- 
cient in continuum descriptions [J, H, H, 0, H, @. The 
potential is truncated at r c = 2.5<r. Newton's equations 
are integrated with a velocity Verlet algorithm[ll| with 
a time step St = 0.005r, where r = erfm/e) 1 / 2 represents 
the characteristic time scale with fluid particle mass m. 
A Dissipative Particle Dynamics (DPD) thermostat [l2T| 
is used to keep the system at a constant temperature 
T = 1.0e/kB, where fee is the Boltzmann constant. A 
coupling constant and a weight function are chosen to 
be C = 0.5r _1 and w(rij) — 1 — respectively, with 
the same cutoff r c = 2.5a as the LJ potential (for fur- 
ther details see Ref.[12]). A Langevin thermostat [l3T] is 
also considered to check the correctness of our results. 
The thermostat with damping constant O.lr -1 is only 
applied in the axial direction of the cylinder in which 
the fluid is not being sheared. We have checked that 
both thermostats give the consistent results confirming 
the appearance of the velocity inversion phenomenon in 
our system. Thus, in this paper, we only present results 
obtained with the DPD thermostat. 

The fluids are confined to a gap between two concen- 
tric cylinders, whose walls are composed of particles of 
mass m w — 100m, tethered by a stiff linear spring with 
constant k w — 100 to fixed lattice sites. The fluid par- 
ticles interact via LJ potential with the Aff = 1. The 
channel wall are made of the same material with the par- 
ticle and interact with fluid particle with Af w = 1 at the 
inner cylinder and Af w = at the outer cylinder. Pe- 
riodic boundary condition is used in the axial direction. 
All simulations are performed with a fixed inner cylinder 
radius as Rija — 2.05 (4 molecular diameters) and dif- 
ferent radius of the outer cylinder ranging from 10.26 (20 
molecular diameters) to 27.36 (53 molecular diameters). 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Density and tangential velocity profiles 

The fluid is sheared by rotating the inner cylinder at a 
constant angular velocity lu = 0.1 rad/r, while the outer 
cylinder remain stationary for all the simulations. After 
an equilibration time of 5 x 10 4 t, the averages are ob- 
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FIG. 2: (Color online) Density profiles for the angular veloc- 
ities ui — O.Orad/r (blue dotted line) and lu — 0.1 rad/r (red 
solid line) with different fluid densities 0.04, 0.40, 0.52, and 
0.82. 



taincd by dividing the interior of the tube into cylindri- 
cal shells of thickness cr/4 for a period of about 2 x 10 5 r. 
Figure [2] shows the density profiles for both the station- 
ary (lu = 0.0 rad/r) and the rotating (lu = 0.1 rad/r) 
inner cylinder denoted by the blue dashed and the red 
solid line, respectively. The fluid densities used in Fig. [3] 
which are volume-average number densities of the fluids, 
p/cr 3 , are 0.04, 0.40, 0.52, and 0.82. There are fluid layers 
near the inner cylinder corresponding to density peaks 
as a consequence of the attractive interaction (Af w = 1) 
between the fluid particles and the inner cylinder, while 
near the outer cylinder with Af w = the peak does not 
exist except for dense regime (see Fig. [5Jd)). For dense 
regime, however, a large number of layers are observed 
near the outer cylinder due to the highly packing of the 
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FIG. 3: Tangential velocities of the flows as a function of the 
radial positions at different volume-average number densities 
of the fluids. The inset shows the slip velocities as a function 
of the fluid density, which are measured at the closest bin to 
the outer cylinder. 
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FIG. 4: Radial force profiles of the flow (circle) for the in- 
dicated values of the fluid density p/a 3 (0.04,0.40,0.52 and 
0.82) as a function of the radial position at uo — 0.1 rad/r. The 
tangential velocity (dotted line) profiles in Fig.[3]are added to 
compare them with the radial force profiles. The horizontal 
solid lines are used as guidelines to the eye. 
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FIG. 5: Tangential velocity profiles with different Af m ranging 
from 0.0 to 0.6. The fluid density is fixed at p/a 3 = 0.40. The 
inset shows the slip velocity near the stationary outer cylinder 
as a function of Af w . 

B. Radial force profile 



particles. When the inner cylinder is rotated, the an- 
gular momentum attained from the rotating cylinder is 
continuously transferred to the outer region through the 
inter-particle interactions, so that the particles near the 
inner cylinder are pushed out to the outer region. 

We have carried out the simulations at different fluid 
densities ranging from 0.04 to 0.82 to measure the tan- 
gential velocity of the flows. The tangential velocity ini- 
tially decreases with the radial position due to the mo- 
mentum dissipation through inter-particle interactions as 
shown in Fig. [3J It, however, is shown that the veloci- 
ties are inverted at below a certain fluid density about 
p/a 3 w 0.6. In other words, the velocity increases with 
distance from the rotating inner cylinder. This unusual 
behavior of the tangential velocity is similar to the results 
of the previous studies [l], 0, H, 0, Q • The velocities finally 
become decaying again near the outer cylinder, which will 
be explained in detail below. It is well known that the 
slip velocity at the stationary outer cylinder plays cru- 
cial role of the velocity inversion [lj. In our simulations 
such slip velocity appear at below a certain fluid den- 
sity p/a 3 « 0.6, below which the velocities are inverted. 
Above 0.6, no velocity slips and no tangential velocity is 
inverted even though the fluid-wall interaction is spec- 
ular (A f w — 0) . Inset of Fig. [3] shows the slip velocity 
near the outer cylinder as a function of the fluid den- 
sity. The maximum value of the slip velocity is located 
at p/a 3 = 0.40, where the degree of the velocity inversion 
is largest. For low density, the slip velocity increases with 
fluid density where the interactions of the particles with 
the wall and with each other are dominant. However, for 
dense cases, it decreases where viscosity becomes domi- 
nant, so that the flow gets slow down. It is sure that the 
velocity inversions are shown if the slip velocity is higher 
than a certain value. 



In order to gain an insight on such an unusual in- 
verted velocity profile in the case of a stationary outer 
cylinder and a rotating inner cylinder, we measure the 
radial force generated by both the rotating inner cylin- 
der with Af w — 1 and the stationary outer cylinder with 
Afw = 0. Figure d] shows the radial force profiles at dif- 
ferent fluid densities, corresponding to 0.04, 0.40, 0.52, 
and 0.82. The tangential velocity profiles in Fig. [3] are 
added to compare them with the radial force profiles. It 
is important to put our finger on the appearance of the 
peaks in the radial force profiles, as in the density profiles 
(see Fig. [5]). However, compared to the peaks of the den- 
sity, the oscillations of the radial force are shown more 
clearly. In the vicinity of the inner cylinder, oscillations 
of the radial force reveal due to the adsorption of the 
fluid at the cylinder. On the other hand, near the reflec- 
tive outer cylinder there exist no such oscillation except 
for the dense regime as in Fig. |3|d). I n Figs. EJa), (b), 
and (c), there is a point at which the outward motion 
of the fluid meets to the inward motion, which results 
from the competition between the centrifugal force by 
the rotating inner cylinder and the centripetal force by 
the reflective outer cylinder. Moreover, due to the lack of 
fluid- wall interaction near the outer cylinder, the angular 
momentum attained from the rotating inner cylinder is 
well transferred to the outer region and induces the ve- 
locity slip. It is interesting to note that the tangential 
velocity begin to be inverted at the point. However, there 
is no such a point for dense case as shown in Fig. H^d) . 
In this regime, the fluid-wall interaction near the outer 
cylinder are very strong, so that the fluids in the first 
layer can hardly slip and the first fluid layer induces a 
second layer by the fluid-fluid interaction. The second 
layer induces a third, and so on. Based on these obser- 
vations, it is sure that the presence of the layering in the 
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FIG. 6: Density (a) and radial force(b) profiles of the flows as 
a function of the radial position at uo — 0.1 rad/r with several 
values of Af w . The fluid density is fixed at p/a 3 = 0.40. The 
dotted line in Fig. (b) denotes a zero value of the radial force 
at which the sign of the radial force is changed. 



vicinity of the outer cylinder is a key factor to determine 
whether the velocity inversion phenomenon appears or 
not. 



C. Layering effect 

To understand the effect of the layering near the outer 
cylinder on the velocity inversion phenomenon, we mea- 
sure the tangential velocity with different values of inter- 
action parameter Af w between the fluid and the outer 
cylinder at fixed fluid density p/cr 3 — 0.4. At this den- 
sity, the velocity inversion behavior can be seen most 
clearly (see Fig. [3]). As Af w increases, the degree of the 
velocity inversion gradually decreases and finally disap- 
pears at exceeding Af w = 0.3 as shown in Fig. \5[ The 
slip velocity shown in the inset of Fig. [5] for six inde- 
pendent realizations of different Af w , corresponding to 
0.0, 0.1, 0.2, 0.3, 0.4, and 0.6 also decreases to zero. These 
indicate that the velocity inversion behavior is strongly 
related to the strength of the fluid-wall interaction near 
the outer cylinder. The responses of the density and the 
radial force with increasing the value of A f w is shown in 
Fig. As the value of Af w increases, the outer cylinder 
becomes more attractive so that more fluid particles can 
move to the outer region and the adsorbed layers begin 
to form near the outer cylinder. This formation of the 
adsorbed layers can be seen most clearly in the radial 
force profile rather than in the density profile as revealed 
in Fig. [6] Even for Af w — 0.0, an indication of a peak of 
the radial force is shown. The fluid particles of the first 
layer almost stick to the wall, which causes the flow to 
retard in the vicinity of the wall. The rapid decay of the 
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FIG. 7: Tangential velocity profiles with four different radius 
of the outer cylinders at given fluid density p/a 3 — 0.40 and 
fixed angular velocity u = 0.1 rad/r. 



tangential velocity near the outer cylinder shown in Fig. [5] 
can be explained by this layering effect even though it is 
weak. When the value of Af w is small (less than 0.3), 
the first layer adsorption is only appeared weakly and its 
size gradually increases as Af w increases from zero. For 
small value of Af w , the velocity of the flow can largely 
slip and the velocity can be inverted. As Af w increases 
beyond 0.3, the effect of the first layer is getting strong 
and induces a second layer, the second layer induces a 
third layer, and so on. From the layer adsorption point 
of view, the fluid particles of a second layer are retarded 
by the first layer, a third layer is retarded by the second 
layer, and so on. As the number of fluid layer increases, 
the no-slip nature at fluid-wall interface spreads into the 
bulk, so that the flow can not slip and no velocity inver- 
sion appears. It is finally concluded that the appearance 
of the velocity inversion phenomenon is determined by 
whether the layer adsorption near the outer cylinder is 
present or not. 

The response of the velocity inversion with increasing 
the radius of the outer cylinder ranging from R/a = 10.26 
to 27.36 at fixed fluid density 0.40 and angular velocity 
0.1 rad/r is shown in Fig. [7] At R/a = 10.26, the tan- 
gential velocity is fully inverted. As the radius of the 
outer cylinder increases, the effect of rotating cylinder is 
reduced, so that finally the velocity would not slip any 
more and not be inverted at given angular velocity. We 
find also that no inverted velocity profile appears in the 
case where the inner cylinder is stationary whereas the 
outer cylinder is rotating. 



IV. CONCLUSIONS 

Using molecular dynamics simulation, we study the be- 
havior of the nanoscale cylindrical Couette flow between 
two concentric rotating cylinders. In simulations, both 
DPD and Langevin thermostats which are suitable for 
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out-of-cquilibrium simulations are used to keep the sys- 
tem at a constant temperature. When the inner cylinder 
rotates and the outer cylinder is at rest, we find that the 
tangential velocity is inverted with a large velocity slip at 
the outer cylinder. The large slip at the outer cylinder, 
which plays an important role in the appearance of the 
inversion, occurs at low fluid density and weak fluid-wall 
interaction near the outer cylinder. We also find that the 
appearance of the inversion behavior is strongly related 
to the presence or absence of the layering of the fluid near 



the outer cylinder. The formation of the layers can be 
seen most clearly in the radial force profile. 
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